Condensation Energy and Spectral Functions in High Temperature Superconductors 
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If high temperature cuprate superconductivity is due to electronic correlations, then the energy difference 
between the normal and superconducting states can be expressed in terms of the occupied part of the single 
particle spectral function. The latter can, in principle, be determined from angle resolved photoemission 
(ARPES) data. As a consequence, the energy gain driving the development of the superconducting state is 
intimately related to the dramatic changes in the photoemission lineshape when going below T c . These points 
are illustrated in the context of the "mode" model used to fit ARPES data in the normal and superconducting 
states, where the question of kinetic energy versus potential energy driven superconductivity is explored in 
detail. We use our findings to comment on the relation of ARPES data to the condensation energy, and to 
various other experimental data. In particular, our results suggest that the nature of the superconducting 
transition is strongly related to how anomalous (non Fermi liquid like) the normal state spectral function 
is, and as such, is dependent upon the doping level. 

PACS numbers: 74.25.-q, 74.25.Bt, 79.60.Bm 



I. INTRODUCTION 

The origin of high temperature superconductivity in 
the cuprates is still a matter of great debate. Recently, 
there have been several different theoretical proposals 
for the mechanism of high T c superconductivity, each of 
which leads to a characteristically different reason for the 
lowering of the free energy. This has focused attention on 
how various spectroscopic probes can yield information 
on the source of the condensation energy which drives 
the formation of the superconducting ground state. 

The first, and perhaps the most radical, proposal is the 
interlayer tunneling theory of Anderson and co-workers, 
where it is conjectured that the condensation energy is 
due to a gain in the c-axis kinetic energy in the supercon- 
ducting state:. Some measurements of the c-axis pene- 
tration depthcl are in conflict with the predictions of this 
theoxjy. Others, such as recent c-axis optical conductivity 
datal3 indicating a violation of the optical sum rule, are 
in support of this hypothesis, although alternative expla- 
nations have been proposed for these observations!! An 
even more unusual suggestion has been recently made by 
Hirsch and Marsiglio,B where they argue that the bulk 
of the condensation energy comes from a gain in the in- 
plane kinetic energy. A rather different approach pro- 
poses the lowering of the Coulomb energy in the long 
wavelength, infrared regionQ, which has not been exper- 
imentally tested as yet. A fourth approach advocates 
a lowering of the exchange energy in the superconduct- 
ing state due to the formatiog.-pf a resonant mode in 
the dynamic spin susceptibilityulS near q = (tt, 7T, 7t), and 
has recently received experimental support from neutron 
scattering studiescl 

We note that all of the above proposals focus on a 
part of the Hamiltonian describing the system: either a 



part of the kinetic energy, or a part of the interaction 
energy. Correspondingly, the experiments to test these 
ideas focus on two-particle correlation functions in a spe- 
cific region of momentum and frequency space. 

In this paper we propose to exploit a very general exact 
relation between the one-particle Green's function of a 
system and its internal energy (see Eq. 1 below) . This ap- 
proach, in principle, allows us to determine the "source" 
of the condensation energy without making any a priori 
assumptions about which piece of the Hamiltonian is re- 
sponsible for the gain in condensation energy. The exact 
expression used involves moments of the occupied part of 
the one-electron spectral function, and since this quan- 
tity is directly relatecO to angle-resolved photoemission 
spectroscopy (ARPES) measurements, our approach also 
appears very promising from a practical point of view. 

As a specific illustration of this general framework, 
we study the condensation energy for a very simple self- 
energy for the normal and superconducting states which 
captures the essential features of the ^observed ARPES 
lineshapes, the so-called mode modeltlrt3. This anal- 
ysis leads to several interesting conclusions as discussed 
below, but most importantly, it suggests an intimate con- 
nection between the nature of the normal state spectral 
function (Fermi liquid or non Fermi liquid) and the for- 
mation of sharply defined quasiparticle excitations below 
T c , and the gain in free energy in the superconducting 
state. 

This paper is organized as follows. In Section II, the 
formalism relating the condensation energy to the spec- 
tral function is developed. In Section III, the mode model 
is introduced, and the nature of the resulting condensa- 
tion energy is discussed. In Section IV, our observations 
concerning ARPES spectra are used to comment on the 
results of previous spectroscopic studies, as well as the 
origin of the condensation energy. In Section V, we ad- 
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dress the question of the nature of the superconducting 
transition versus hole doping. In Section VI, we offer 
some concluding remarks. Finally, we include two ap- 
pendices. Appendix A further explores questions raised 
in Section II in regards to the full Hamiltonian and the 
virial theorem. In Appendix B, we comment on the ap- 
plicability of the formalism of Section II to experimental 
data (ARPES and tunneling). 

II. FORMALISM 

We begin with the assumption that the condensa- 
tion energy does not have a component due to phonons, 
though as we mention below, this condition can be re- 
laxed. We note that at optimal doping, the isotope ex- 
ponent a is essentially zercO, and Chesteillj proved that 
the change in ion kinetic energy between superconduct- 
ing and normal states vanishes for a = 0. To proceed, 
we assume an effective single-band Hamiltonian which 
involves only two particle interactions. Then, simply ex- 
ploiting standard formulasE§liII for the internal energy 
U = (H — /iN) ([i is the chemical potential, and N the 
number of particles) in terms of the one-particle Green's 
function, we obtain 

U N -U S = 

/+oo 
du(u + efc)/H [Ajv(k, La) - A s (k, lj)] (1) 

Here and below the subscript N stands for the normal 
state, S for the superconducting state. j4(k, cj) is the 
single-particle spectral function, f(oj) the Fermi function, 
and efc the bare energy dispersion which defines the ki- 
netic energy part of the Hamiltonian. Note that the fj,N 
term has been absorbed into u> and e^, that is, these 
quantities are defined relative to the appropriate chemi- 
cal potential, /ijy or ^5. In general, /ijy and us will be 
different. This difference has to be taken into account, 
since the condensation energy is small. 

The condensation energy is defined by the zero tem- 
perature limit of Un — Us in the above expression. Note 
that this involves defining (or somehow extrapolating to) 
the normal state spectral function at T = 0. Such an ex- 
trapolation, which we return to below, is not specific to 
our approach, but required in all estimates of the conden- 
sation energy. We remark that Eq. 1 yields the correct 
condensation eneijgy, iV(0)A 2 /2, for the BCS theory of 
superconductivity!!^. 

We also note that Eq. 1 can also be broken up into 
two pieces to individually yield the thermal expectation 
value of the kinetic energy (using 2efc in the parentheses 
in front of /(w)), and that of the potential energy (using 
to — ek instead). Further, this expression can also be 
generalized to the free energy by including the entropy 
term as discussed by Wadalia. Moreover, if the phonons 
can be treated in an harmonic approximation, the terms 



missing in Eq. 1 (half the electron-phonon interaction, 
and all other phonon terms) reduce to twice the phonon 
kinetic energy! 1 J£3. The phonon kinetic energy can. then 
be determined if the isotope coefficient is knownEj. For 
a = 1/2, the missing terms in this approximation reduce 
to twice the condensation energy, so that Eq. 1 is realized 
again, but with a negative sign. 

The great advantage of Eq. 1 is that it involves just 
the occupied part of the single particle spectral function, 
which is measured Jay angle resolved photoemission spec- 
troscopy (ARPES)E9. Therefore, in principle, one should 
be able to derive the condensation energy from such data, 
if an appropriate extrapolation of the normal state spec- 
tral function to T=0 can be made. On the other hand, 
a disadvantage is that the bare energies, e^, are a pri- 
ori unknown. Note that these are not directly obtained 
from the measured ARPES dispersion, which already in- 
cludes many-body renormalizations, nor are they sim- 
ply determined by the eigenvalues of a band calculation, 
as such calculations also include an effective potential 
term. Rather, they could be determined by projecting 
the kinetic energy operator onto the single-band sub- 
space. Methodologies for doing this when reducing to an 
effective single-bandJIubbard model have been worked 
out for the cupratescJ, and could be exploited for this 
purpose. 

Eq. 1 trivially reduces to the following: 
U N - U s = £ fe M k ) ~ n s (k)] 

k 

/+00 
dww/(w) [N N (u) - N s (lj)} (2) 
-00 

where n(k) is the momentum distribution function, and 
N(ui) the single-particle density of states. While ARPES 
has the advantage of giving information on both terms 
in this expression, other techniques could be exploited as 
well for the individual terms in Eq. 2. For instance, n(k) 
in principle can be obtained from positron annihilation 
or Compton scattering, while N(u>) could be determined 
from tunneling data, although matrix elements could be 
a major complication for both tunneling and ARPES. 

We conclude this Section with some remarks about 
a low-energy effective single-band Hamiltonian used to 
derive Eq. 1 versus the full Hamiltonian of the solid 
which includes quadratic dispersions for all (valence and 
core) electrons, ionic kinetic energies, together with all 
Coulombic interactions (see e.g., Ref. As shown by 
ChesteitU the full H can be very useful for studying the 
condensation energy. We discuss some points related to 
such a description in the Appendix. 

Here we only wish to emphasize one important point 
which will come up later in our analysis. In terms of the 
full Hamiltonian, the transition to the superconducting 
state must be driven by a gain in the potential energy 
(ignoring ion terms for this argument), as is intuitively 
obvious and also rigorously shown by Chester using the 
virial theorem. However, the kinetic energy terms in 
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the effective single-band Hamiltonian can (and in gen- 
eral do) incorporate effects of the potential energy terms 
of the full Hamiltonian. Further, there is no virial the- 
orem restriction on the expectation values of the kinetic 
and potential terms of the effective Hamiltonian (since 
these do not, in general, obey the requisite homogeneity 
conditions). As a consequence, there is nothing prevent- 
ing the effective low-energy Hamiltonian from having a 
superconducting transition driven by a lowering of the 
(effective) kinetic energy. 



III. THE MODE MODEL 

To illustrate the power of the formalism, as well as 
some of the subtleties discussed above, we now analyze 
the condensation energy arising from a spectral function 
described by a simple model self-energy which captures 
some of the essential features of the ARPES data in the 
important region of the Brillouin zone near (w, 0) in the 
cuprates. These features are: (1) a broad normal state 
spectral function A which seems T-independent in the 
normal state (except in the underdoped case, where there 
is a pseudogap which fills in as T increases) , and thus can 
be used as the extrapolated "normal" state An down to 
T = in Eq. 1. (2) A superconducting state spectral 
function As which shows a gap, a sharp quasiparticlc 
peak, and a dip-hump structure at higher energies. At 
a later stage, we will have to make some reasonable as- 
sumptions about the k-dependence of the spectral func- 
tions to perform the zone sum in Eq. 1. 

These nontrivial changes in the ARPES lineshape go- 
ing from the normal to the superconducting state have 
been attributedEilllj to the interaction of an electron with 
an electronic resonant mode below T c , which itself arises 
self-consistently from the lineshape change. Strong argu- 
ments have been given which identify this resonant mo Hp 
with one observed by magnetic neutron scatteringtilEj. 
Thus our analysis below will also have bearing upon the 
arguments mentioned in the Introduction which relate 
the resonant mode directly to changes in the exchange 
energy. 

The simplest version of the resonant mode model is a 
self-energy of the form 



E = -In 

7T 



UJ — UJq — A 



- uj + A 



+ »re(M - a* - A), (3) 



where ljq is the resonant mode energy, A the supercon- 
ducting energy gap, and the step function. (A more 
complicated form has been presented in earlier workE3.) 
This self-energy is-then used in the superconducting state 
spectral functiontj 



A 1 T ^ W 

A = - lm ^2T^2 



(4) 



where Z = 1 — E/w. We note that for this form of E, the 
spectral function As will consist of two delta functions 



located at ±E, where E satisfies two conditions: (1) it 
has a value less than lj q + A and (2) the denominator of 
Eq. 4 vanisheSp-The weight of the delta functions are then 
determined asE3 \dA~ 1 {±E) j 'du>\. In addition, there are 
incoherent pieces for \oj\ greater than loq + A. We use the 
same self-energy for the (extrapolated) normal state with 
A = and luq = 0, so that Am reduces to a Lorentzian 
centered at e with a full width half maximum of 2I\ 
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FIG. 1. (a) Condensation energy contribution, E con d, ver- 
sus single particle energy e for the model self-energy of Eq. 3. 
As discussed in the text, the quantity plotted is the result (as 
a function of e) after the ^-integration is done in Eq. 1. The 
parameters are T=230 meV, A=32 meV, and Wo=41.6 maW, 
which were obtained from fits to ARPES spectra at (ty, 0)l3. 
The normal state is obtained by setting Uq and A to zero. 
The dotted lines are a decomposition of E con d into separate 
kinetic and potential energy pieces, (b) Condensation energy 
contribution for the BCS theory using the same A. (c) A 
repeat of (a), but with the superconducting state replaced by 
the normal state with wo=41.6 meV (and so labeled as A_E 
instead of E cond ). 
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To begin with, for simplicity, we treat both loq and 
A as momentum independent. It is straightforward to 
evaluate Eq. 1 with the sum over momentum reducing to 
an integral over e. In Fig. la, we plot the integrand of the 
e integral (i.e., after the u> integral has been doricjl. The 
parameters used are the same ones used earlienij to fit 
ARPES data near optimal doping at the (ir, 0) point. The 
result is somewhat surprising. The integrand is negative 
for e near zero (i.e., k near kp) and positive for e far 
enough, away. This should be contrasted with the BCS 
resultfij, shown in Fig. lb, where the contribution at kp 
(which is A/2) is maximal and positive. 

To gain insight into this unusual result, we also show in 
Fig. la the decomposition of this result into kinetic and 
potential energy pieces. Unlike BCS theory (Fig. lb), 
where the condensation is driven by the potential en- 
ergy, in the mode model case, it is kinetic energy driven. 
To understand the unusual decrease in the kinetic energy 
as one goes below T c , we show in Fig. 2 the momentum 
distribution function n(k) plotted versus e. Note that 
in contrast to BCS theory, n(k) is sharper in the super- 
conducting state than in the normal state. The reason 
is very simple. The (extrapolated) normal state is sub- 
ject to a large broadening T all the way down to T = 
which smears out n(k) on the scale of T. At T = the 
result is simply: n.Ar(k) = 1/2 — tan -1 (e/T) /ir. In the 
superconducting state, although A broadens n(k) as in 
BCS theory, one now has quasiparticle peaks. The ef- 
fect of this on sharpening n{k) is much larger than the 
broadening due to A (for A <C T), so the net effect is 
a significant sharpening. As a consequence, the kinetic 
energy is lowered in the superconducting state. 
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FIG. 2. Momentum distribution function versus e in the 
superconducting state (SC), normal state with cjo=0 (NS), 
and in the normal state with cjo=41.6 meV (NS-mode). Same 
parameters as Fig. 1. 



Note that these counterintuitive results would not have 
been obtained had loq retained the same (non-zero) value 
in the normal state. In this case, sharp quasiparticles 
would exist in the normal state, and all of our usual ex- 
pectations are fulfilled: njv(k) would have had a step 
discontinuity (also illustrated in Fig. 2), and the normal 
state kinetic energy would have been considerably lower 
than the superconducting one. In fact, for this situation, 
the model is equivalent to that of Einstein phonons in an 
approximation where the gap is treated as a (real) con- 
stant in frequency^. However, the normal state ARPES 
data near (jr, 0) are clearly consistent with uq = and 
are T-independent with a T » T, which suggests that 
the T=0 extrapolation used here is reasonable. 

These points are further illustrated in Fig. lc, where 
we show the energy difference between the normal state 
with loq non-zero, and the normal state with loq zero. 
Note the similarity to Fig. la, i. e. , the unusual behavior 
in Fig. la is due to the formation of a gap in the inco- 
herent part of the spectral function, with the resulting 
appearance of quasiparticle states, and thus not simply 
due to the presence of a superconducting energy gap, A. 
Now, in the real system, it is the transition to a phase 
coherent superconducting state which leads to the ap- 
pearance of the resonant mode at non-zero energy, which 
causes the gap in 7m£, which results in the incoherent 
gap and quasiparticles, which in turn generates the mode. 
Although this self-consistency loop clearly indicates the 
electron-electron nature of the interaction (as opposed 
to an electron-phonon one), the connection of these ef- 
fects with the onset of phase coherence (as opposed to 
the opening of a spectral gap, which is known to occur 
at a higher temperature, T*) is not understood at this 
time. That is, the mode model is a crude simulation of 
the consequences of some underlying microscopic theory 
which has yet to be developed. 

As for the potential energy piece, we note that the 
contribution to Eq. 1 at fcp (where e/- = 0) reduces to 
the first moment of the spectral function. In Fig. 3a, we 
plot the spectral function at kp in both the normal and 
superconducting states. (For illustrative purposes, we 
have replaced the delta function peaks in the supercon- 
ducting state by Lorentzians of half width half maximum 
10 meV). From this plot, we note that the quasiparticle 
peaks give a positive contribution to the condensation 
energy, but that at higher energies (large \oj\), there is a 
negative contribution. This negative contribution is very 
important because it is weighted by uj in the integrand 
of Eq. 1. To see this quantitatively, we plot in Fig. 3b 
the first moment difference at the Fermi surface (e& = 0) 
as a function of the lower cut-off on the u> integration 
(the upper cutoff at T = is u = 0). We clearly see the 
positive contribution due to the quasiparticle peak and 
the (five times larger) negative contribution due to the 
incoherent tail. This explains why the net contribution 
from the potential energy term is negative. We can con- 
trast this with BCS theory, where only the quasiparticle 
part exists, and so the net contribution is positive. 
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FIG. 3. (a) Spectral function at the Fermi surface (e = 0) 
in the superconducting (SC) and normal states (NS). (b) First 
moment contribution of (a) to the condensation energy ver- 
sus the lower cut-off in the id integration in Eq. 1. Note posi- 
tive contribution of the quasiparticle peak, and large negative 
contribution from the high energy tail. Same parameters as 
Fig. 1. 

An interesting question concerns what happens in this 
model as the broadening, T, is reduced. In Fig. 4, we 
show results like for Fig. la, but for various V values. As 
r is reduced and becomes comparable to A, one crosses 
over from the unusual behavior in Fig. la to a behavior 
very similar to that of BCS theory in Fig. lb. That is, 
the condensation energy crosses over from being kinetic 
energy driven to being potential energy driven. This is 
not a surprise, since in the limit T goes to zero, the model 
reduces to BCS theory. The physics behind this, though, 
is quite interesting. For large T, the normal state is very 
non Fermi liquid like. As T, is reduced, tbxpgh, the nor- 
mal state becomes more Fermi liquid likecfl. As a conse- 
quence, one crosses over from being kinetic energy driven 
to potential energy driven (when r ~ A). The relation 
of kinetic energy driven behavior with the presence of a 
non Fermi liquid normal state, and a Fermi liquid suaei;- 
conducting state, was realized early on by AndersonOEJ, 
and will be returned to again in Section IV of the pa- 
per. Fig. 4 also draws attention to the fact that being 
kinetic or potential energy driven is a relative point. Note 
in Fig. 4b that near /cp, the two contributions have the 
same sign. Individual terms, such as the potential energy 
in Fig. 4b, and the kinetic energy in other cases we have 
explored, can even change sign as a function of e^. 
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FIG. 4. Condensation energy contribution, as in Fig. la, 
for (a) r=16 meV, (b) 32 meV, and (c) 64 meV. A and uo are 
both 32 meV. Note the crossover from kinetic energy driven 
behavior to potential energy driven behavior as F is reduced. 

A potential worry which emerges from the above cal- 
culations is the large contribution in Figs. 1 and 4 at 
large |e|. In particular, in most cases, the bulk of the 
contribution to the condensation energy comes well away 
from the Fermi surface, in contrast to BCS theory. In 
Fig. la, this is due to the large T, which leads to a sub- 
stantial rearrangement of the spectral function even for 
large |e|, causing large contributions to both the potential 
and kinetic energy pieces. Even in the case of Fig. 4a, 
where T is quite small, there is still a potential energy 
contribution at large |e|. This can be traced to the gap 
in the incoherent part of the spectral function, with the 
resulting spectral weight being recovered around lo = e, 
leading to a potential energy shift. Even in the BCS 
case, Fig. lb, the individual potential and kinetic energy 
pieces would not converge if integrated over an infinite 
range in e (even though their difference would). In BCS 
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theory, this is corrected by an ultraviolet cut-off (the De- 
bye frequency). We elect not to include such a cut-off in 
the mode model, since it would be lead to another ad- 
justable parameter, and the e integral is bound by the 
band edges, and so is convergent. In the real system, 
the "mode" effects in the spectral function disappear as 
one approaches the band edges, and as discussed in the 
following paragraph, this effect can be crudely simulated 
by setting the mode energy proportional to A/., the lat- 
ter quantity in the d-wave case vanishing along the zone 
diagonal where the band edges are located. 

Although we plot only the differences in Figs. 1, 3b, 
and 4, the individual normal and superconducting state 
terms are quite large. This raises the question of what the 
value of Eq. 1 would actually be if summed over the zone. 
To do this, we must make some assumptions about what 
the momentum dependence of various quantities are. For 
simplicity sake, we will treat r as k independent, though 
we note that available ARPES data are consistent with 
this quantity being reduced in size as one moves from 
(-7T, 0) towards the Fermi crossing along the (tt, tt) direc- 
tion. In the first sum, denoted by case (a), we treat A and 
luq as k independent. In the second sum, denoted by case 
(b), we replace A by A*, = Aofcos^a) — cos(fc y a))/2, 
where A& is the standard d-wave gap function, but still 
retain a fc-independent loq- In the third sum (c), in addi- 
tion to the d-wave A*, we also take ujq — c|Afc|, with the 
k dependence of ojq crudely simulating the fact that the 
mode effects in the spectral function are reduced as one 
moves away from the (tt, 0) points of the zoneEl The val- 
ues of these parameters are the same as used in Fig. la, 
and are consistent with ARPES and neutron data for 
Bi2212 (r=230 meV, A =32 mcV, c=1.3). To perform 
the zone sum, we have to make some assumptions on 
what the ej, are. As the mode model is designed to ac- 
count for the difference between the normal state and su- 
perconducting state, we elect to use normal state ARPES 
dispersions for ejJHlJ, though we caution that this repre- 
sents a different choice for the "kinetic" energy part of 
the effective single-band Hamiltonian than is typically 
usedcJ. Because this dispersion has particle-hole asym- 
metry, the chemical potential will not be the same in the 
superconducting state as in the normal state. The chem- 
ical potential is thus tuned to achieve the same density (a 
hole doping x=0.16) as the normal state. Note that the 
normal state density itself is a function of T (we assume 
u}o=0 for the normal state). 

Performing the zone sum, we find condensation ener- 
gies of +3.6, +3.3, and +1.1 meV, per CuO plane, for 
case (a), (b) and (c) respectively. We note that the last 
result is the more physically appropriate, and though 
small, is somewhat larger than the condensation energy 
of 1/4 meV per plane estimated by Loraroje/ al from spe- 
cific heat data for optimal doped YBCOEj. The above 
values will be reduced if a more realistic k, ui dependence 
is used for T, since, as we noted above, V decreases as 
one moves away from (tt, 0). As consistent with Fig. la, 
the contribution to the condensation energy is negative 



for an anisotropic shell around the Fermi surface (due to 
the anisotropy of A& and e/c), and positive outside of this 
shell. Again, this will be sensitive to the k dependence 
of r, as can be seen from Fig. 4. We also remark that 
there are chemical potential shifts of +2.6, +2.1, and 
+1.4 meV, respectively, for case (a), (b) and (c). Again, 
the last value is the more physically appropriate. It is 
very interesting to note that somewhat smaller positive 
shifts (around +0.6 meV) have been seen experimentally 
in YBCOEj. These shifts are a consequence of particle- 
hole asymmetry and the change in rik when going into 
the superconducting state. 

IV. CONNECTIONS WITH PREVIOUS WORK 

While a quantitative evaluation of Eq. 1 using experi- 
mental data as input on the right hand side must await 
further progress as discussed in Appendix B, several qual- 
itative points can be made even at this stage. From 
Eq. 1, there is a one to one correspondence between the 
changes in the spectral function and the condensation 
energy. That is, the condensation energy is due to the 
profound change in lineshape seen in photoemission data 
when going below T c . When summed over the zone, this 
in turn leads to changes in the tunneling density of states 
(second part of Eq. 2). These spectral function changes 
cause, and are themselves caused by, changes of vari- 
ous two particle correlation functions, such as the optical 
conductivity and the dynamic spin susceptibility, which 
have previously been used by others to comment about 
the nature of the condensation energy. 

In this context, we now discuss the earlier work con- 
cerning the c-axis conductivity. The most dramatic 
changes in the ARPES lineshape when going below T c 
occur near the (tt, 0) points of the zone. It is exactly 
these points of the zone which appear to have the largest 
c-axis tunneling matrix elements associated with themEil. 
Previous work has found a strong correlation between the 
c-axis conductivitjj— and ARPES spectra near the (tt, 0) 
points of the zonou. Therefore, it is rather straight- 
forward to speculate that it is the formation of strong 
quasiparticle peaks in these regions of the zone, and the 
resulting changes in the spectral function at higher bind- 
ing energy, which is responsible for the lowering ofJici 
c-axis kinetic energy. We note that earlier, AndersonEj-El 
had remarked that if the quasiparticle weight is coming 
from high binding energy, then one would expect a lower- 
ing of the kinetic energy. This in fact is what is occuring 
in the mode model calculations, though we note from our 
work that the true quantity which determines the sign of 
the kinetic energy change in the vicinity of kp is the gra- 
dient of the momentum distribution function at kp. 

We also remark that the change in c-axis kinetic en- 
ergy has been recently addressed by Ioffe and Millis in 
the context of the same mode model used in the current 
papercl. These effects would enter directly in Eq. 1 by 
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including a c-axis tunneling contribution to As for 
the in-plane kinetic energy, it is so large that it is difficult 
to determine its contribution to Eq. 1 from optical con- 
ductivity data because of some of the same normalization 
concerns mentioned in Appendix B in regards to ARPES 
and tunneling data. Still, if the mode model calculation 
is a reflection of reality, we can speculate that n(k) will 
probably sharpen in the superconducting state, leading 
to a lowering of the in-plane kinetic energy. How large 
the effect will be is somewhat difficult to determine, in 
that the same regions of the zone where large changes are 
seen in the ARPES lineshape are also characterized by 
small Fermi velocities (the optical conductivity involves a 
zone sum weighted by vp). Along the (ir, it) direction, for 
instance, there is still some controversy concerning how 
dramatic the lineshape change is below r c Eo. Also, 
as can be seen from Fig. 4, this question is very depen- 
dent on the variation of the normal state lineshape in 
the zone. Although the lineshape near (tt, 0) is highly 
non Fermi liquid like, the behavior along the C^rZt) di- 
rection appears to be marginal Fermi liquid likeLjO. As 
remarked in Section III, the more Fermi liquid like the 
normal state lineshape is, the greater the tendency is to 
switch over to potential energy driven behavior instead. 
Improved experimentation should again lead to a resolu- 
tion of these issues. 

This brings us to the question concerning the relation 
of the magnetic resonant mode observed by neutron scat- 
tering to the condensation energy. All calculations of 
the resonant mode assume the existence of quasiparti- 
cle peaks. In the absence of such quasiparticle peaks, 
a sharp resonance is not expected. That is, the sharp 
resonance observed by neutron scattering, and the re- 
sulting lowering in the exchange energy part of the t-J 
Hamiltonian, is again a consequence of the formation of 
quasiparticle states. In this context, it is important to 
note that the d-wave coherence factors associated with 
quasiparticle states are important for the formation of 
the resonance, whether iu the context of calculations in 
the particle-hole channeEj, or in the particle-particle sce- 
nario proposed by Demler and Zhango. In any case, this 
again supports our statement, motivated by Eq. 1, that 
it is the dramatic change in the ARPES spectra below 
T c which is the source of the condensation energy. 

In this regard, we note a puzzling feature in connec- 
tion with the mode model. Although it was designed to 
take into account the effect of the magnetic resonance 
mode on the spectral function, the condensation in the 
mode model is kinetic energy driven. This is in contrast 
to the potential energy driven nature of the condensation 
with the-resonant mode discussed in the context of the t-J 
modelTa, despite the same underlying physics. There are 
two possibilities for this apparent discrepancy. First, the 
break-up of the Hamiltonian into potential and kinetic 
energy pieces depends on the particular single-band re- 
duction which is done. The superexchange energy, which 
is a kinetic energy effect at the level of the Hubbard 
modeEa, appears as a potential energy term when re- 



duced to the t-J Hamiltonian. In the mode model, the 
kinetic energy is equated to based on normal state 
ARPES dispersionsc3, while the potential energy term 
leads to effects described by the E of Eq. 3. 

The second possibility is that the argument of Ref. || 
is confined to low energies of order A. As demonstrated 
in Fig. 3b, if the mode model is confined to such energy 
scales, the first moment (i.e., the potential) term would 
reverse sign, since the quasiparticle peak always gives a 
positive contribution to the first moment. That is, one 
would expect the resonance to lower the exchange energy 
since it is a consequence of the quasiparticle states, which 
lower the potential energy in Eq. 1. It is the difference in 
the high energy incoherent tails (Fig. 3), though, which is 
ultimately responsible for the increase of the net potential 
energy in Fig. la. This would imply that the neutron 
scattering resultsu may change if more complete data at 
higher energies and other q values are obtained. That is, 
the true answer will depend on where the weight for the 
neutron resonance is coming from, in complete. analogy to 
the earlier mentioned argument of AndersorcJ in regards 
to where the quasiparticle weight is coming from. 

This discussion again emphasizes that the current de- 
bate concerning kinetic energy driven superconductivity 
versus potential energy driven superconductivity must be 
kept in proper context, as the very definition of the ki- 
netic and potential pieces is dependent upon what ef- 
fective low energy Hamiltonian one employs, and what 
energy range one considers. 



V. DOPING DEPENDENCE 

The condensation energy as estimated from specific 
heat isr-known to decrease strongly as the doping is 
rcducjcdp. This is despite the increase of the spectral 
gapEIlHEj. There are two reasons for this suggested by 



the above line of reasoning. First, the normal state itself 
at T c already exhibits a large spectral gap, the so-called 
pseudogap, which acts to reduce the difference in Eq. 1. 
Second, the weight of the quasipactiple peak strongly de- 
creases as the doping is reducedli3. This reduces the 
quasiparticle contribution to both the first moment and 
to n(k). We caution that the normal state extrapola- 
tion down to T=0 will be more difficult to estimate for 
underdoped experimental data because of the influence 
of the pseudogap, which is known to fill in as a func- 
tion of temperatureO. Still, the available underdoped 
ARPES and tunneling data are certainly in support of 
a smaller condensation energy than overdoped data due 
to the pseudogap, which is in agreement with conclu- 
sions based on specific heat dataO. The new contribu- 
tion to these arguments is the strong reduction of the 
weight of the quasiparticle peak in the underdoped case 
which makes the condensation energy smaller still. In 
fact, based on our arguments, the strong reduction of the 
superfluid density upon underdoping is almost certainly 
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connected with the strong reduction in the quasiparticle 
weight. j— - 

Finally, Andersoncil has speculated that the supercon- 
ducting transition temperature is potential energy driven 
on the overdoped side, kinetic energy driven on the un- 
derdoped side. This is a distinct possibility, since T is 
known from ARPES dataE3 to be strongly reduced as 
the doping increases on the overdoped side, and as Fig. 4 
demonstrates, one might expect (if the mode model is 
a reflection of reality) a crossover from kinetic energy 
driven behavior to potential energy driven behavior as U 
is reduced. In this context, we note the Basov et aZresultcl 
that the lowering of the c-axis kinetic energy appears to 
be confined to the underdoped side of the phase diagram. 
Moreover, if one attributes T* oit-fche underdoped side to 
the onset of pairing correlational, then one anticipates 
a potential energy gain due to pairing to occur at this 
finite temperature crossover. At T c , phase coherence in 
the pair field is-established, and the resulting quasipar- 
ticle formations and related spectral changes could lead 
to a kinetic energy driven transition of the sort discussed 
above. We emphasize "could", since in the context of 
Eq. 1, there is no unambiguous evidence yet from real 
ARPES data that such is the case. 



VI. CONCLUDING REMARKS 

We conclude this paper by noting that the above argu- 
ments based on condensation energy considerations high- 
lights one of the key question of the high T c problem: 
why do quasiparticle peaks only appear below T c ? This 
is especially relevant in the underdoped case, since the 
spectral gap turns on at a considerably higher temper- 
ature thaa— jT c , but the quasiparticle peaks again form 
only at Tji3. This implies that there is a deep connec- 
tion between the onset of phase coherence in the pair 
field and the onset of coherence in the single electron de- 
grees of freedomEl. We suggest that the understanding 
of this connection will be central to solving the high T c 
problem. The result of the current paper is that Eq. 1 
brings this issue into much shaper focus. In particular, 
as a cautionary note, the incoherent part of the spectral 
function is likely to be as important as the quasiparti- 
cle component in determining the condensation energy 
(Fig. 3). That is, it is the overall shape of the spectral 
function (the peak-dip- hump behavior of Fig. 3a), rather 
than just the quasiparticle part, which is ultimately re- 
sponsible for the total condensation energy. We believe 
that experimental data analyzed in the context of Eq. 1 
will play an important role in providing a solution to the 
high T c problem. 
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APPENDIX A: THE FULL HAMILTONIAN AND 
THE VIRIAL THEOREM 

In this Appendix, we make further comments on some 
issues which were briefly discussed at the end of Section 
II, relating to the use of the full Hamiltonian versus an 
effective single-band Hamiltonian. 

We note that as written, Eq. 1 does not apply to the 
full Hamiltonian of the solid which includes all the elec- 
tronic and ionic degrees of freedom together with their 
Coulombic interactions as discussed in Ref. [l5|. In princi- 
ple an expression similar to Eq. 1 could be written if the 
quantities in Eq_l were replaced by matrices in recipro- 
cal lattice spaceuj. For our purposes, where an energy 
difference is being looked at, a unitary transformation 
to band index space would be desirable. The resulting 
off-diagonal terms would then represent interband tran- 
sitions. These could be of potential importance, even for 
the energy difference. For example, the violation of the 
c-axis optical conductivity sum ruletl implies a change 
in interband terms so that the total optical sum rule is 
satisfied. 

The usefulness of t]p,,£ull Hamiltonian is that one can 
use the virial theoremllalB 2K — nV — 3PQ = 0, exploit- 
ing the fact that the kinetic energy K is a homogeneous 
function of order 2 in momentum, and the potential en- 
ergy V is a homogeneous function of order n in position. 
Here P is the pressure and Q, denotes the volume. For 
Coulomb forces n = — 1, and ignoring the pressure terms 
(which are negligible at ambient pressure), this reduces 
to 2K + V = 0. 

If we assume that the form of Eq. I applies to the 
full Hamiltonian (which could be possible if all interband 
terms dropped out of the energy difference, as well as all 
electron-ion and ion-ion terms) then by using the virial 
theorem, the right hand side of Eq. 2 can be shown to 
reduce to 2/3 the first moment of the density of states 
at T = 0. In addition, the change in the kinetic energy 
would be the negative of the condensation energy, with 
the potential energy twice the condensation energy. 

This reduced form of Eq. 2, though, must be treated 
with extreme caution, and is likely not useful to the prob- 
lem at hand. The reason is that the kinetic energy and 



8 



potential energy terms of the full Hamiltonian are not 
the same as the kinetic and potential energy terms of 
the effective single-band Hamiltonian. It is only for the 
former that the virial theorem manipulations would be 
allowed. As an example, BCS theory obeys Eq. 2, but 
not the reduced form. 



APPENDIX B: COMMENTS ON ARPES AND 
TUNNELING 

The purpose of Section III was to demonstrate how 
Eq. 1 works out in practice for a model where exact cal- 
culations could be done. This is important when con- 
sidering real experimental data. We have spent consid- 
erable effort analyzing Eqs. 1 and 2 using experimental 
data from ARPES and tunneling as input, and plan to 
report on these endeavors in a future publication. But 
given what we have learned from the mode model, some 
of the problems associated with an analysis based on ex- 
perimental data can be appreciated. First, the condensa- 
tion energy is obtained by subtracting two large numbers. 
Therefore, normalization of the data becomes a central 
concern. Problems in this regard when considering n(k), 
which is the zeroth moment of the ARPES-data, were 
discussed in a previous experimental paperCJ. For the 
first moment, these problems are further amplified due 
to the cu weighting in the integrand. This can be appre- 
ciated from Fig. 3, where the bulk of the contribution in 
the mode model comes from the mismatch in the high en- 
ergy tails of the normal state and superconducting state 
spectral functions. When analyzing real data, we have 
found that the tail contribution, either from ARPES or 
from tunneling, is very sensitive to how the data are nor- 
malized. Different choices of normalization can even lead 
to changes in sign of the first moment. 

Another concern concerns the k sum in Eq. 1. Both 
ARPES and tunneling have (their own distinct) k- 
dependent matrix elements, which lead to weighting fac- 
tors not present in Eq. 1. For ARPES, these effects can in 
principle be factored out by either theoretical estimates 
of the matrix elementsCj, or by comparing data at dif- 
ferent photon energies to obtain information on themED. 
For tunneling, information on matrix elements can be ob- 
tained by comparing different types of tunneling (STM, 
tunnel junction, point contact), or by employing direc- 
tional tunneling methods. 

Another issue in connection with experimental data is 
an appropriate extrapolation of the normal state to zero 
temperature. Information on this can be obtained by an- 
alyzing the temperature dependence of the normal state 
data, remembering that the Fermi function will cause 
a temperature dependence of the data which should be 
factored out before attempting the T = extrapolation. 
We finally note that the temperature dependence issue is 
strongly coupled to the normalization problem mentioned 
above. In ARPES, the absolute intensity can change due 



sorbed gasses, 
13. In tunnel- 



to temperature dependent changes in at 
surface doping level, and sample locatior 
ing, the absolute conductance can change due to temper- 
ature dependent changes in junction characteristics. In 
both cases, changes of background emission with temper- 
ature is another potential problem. 

Despite these concerns, we believe that with careful 
experimentation, many of these difficulties can be over- 
come, and even if an exact determination of Eq. 1 is not 
possible, insights into the origin of the condensation en- 
ergy will certainly be forthcoming from the data. This is 
particularly true for ARPES, which has the advantage of 
being k resolved and thus giving one information on the 
relative contribution of different k vectors to the conden- 
sation energy. 
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